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Turbulent shear flows, such as those occurring in the wall region of turbulent boundary 
layers, manifest a substantial increase of intermittency with respect to isotropic con- 
ditions. This suggests a close link between anisotropy and intermittency. However, a 
rigorous statistical description of anisotropic flows is, in most cases, hampered by the 
inhomogeneity of the field. This difficulty is absent for the homogeneous shear flow. For 
this flow the scale by scale budget is discussed here by using the appropriate form of the 
Karman-Howarth equation, to determine the range of scales where the shear is dominant. 
The issuing generalization of the four-fifths law is then used as the guideline to extend to 
shear dominated flows the Kolmogorov-Obhukhov theory of intermittency. The procedure 
leads naturally to the formulation of generalized structure functions and the description 
of intermittency thus obtained reduces to the K62 theory for vanishing shear. Also here 
the intermittency corrections to the scaling exponents are carried by the moments of the 
coarse grained energy dissipation fleld. Numerical experiments give indications that the 
dissipation fleld is statistically unaffected by the shear, thereby supporting the conjec- 
ture that the intermittency corrections are universal. This observation together with the 
present reformulation of the theory gives reason for the increased intermittency observed 
in the classical longitudinal velocity increments. 



1. Introduction 

At large Reynolds number, turbulent flows exhibit strong intermittent effects which 
can be described in a number of ways, such as intense coherent structures, anomalous 
scaling in the inertial range, exponential tails in the probability density functions (see 



e.g. Friscb). In the case of homogeneous and isotropic turbulence, more than twenty years 
of intense scientiflc research provide us with strong indications that intermittency shows 
universal properties ( ^reenivasan fc Antonia ). 



As soon as we consider non homogeneous non isotropic turbulence, a well established 
conclusion concerning universal properties of intermittency is still missing, although some 



preliminary results have been proposed. See for instance the papers of Toschi et al. 



Danaila et al., Arad et al.. A particular case of non isotropic turbulence is provided 



by homogeneous shear flows which have been the subject of a number of numerical 



(see Pumir and references therein) and experimental investigations (Shen & War haft) 



Homogeneous shear flows are presumably the simplest non trivial example of non isotropic 
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turbulence where issues such as the generation and the dynamics of coherent vortical 
structures, the effect of the large scale non isotropic forcing on the small scale intermittent 
fluctuations of velocity and the universal (if any) scaling properties of velocity increments 
in the inertial range can be addressed. 

In the present paper we investigate intermittency effects in homogeneous shear flows 



by using a long time highly resolved Direct Numerical Simulation (DNS), see Gualtieri 



et al. for details. Our basic idea is to study the fluctuations of the velocity field by using 
a generalization of the Kolmogorov four-fifths equation. Although for the homogeneous 
shear flow no use of the isotropic constraints can be made, we are still able to write the 



appropriate form of this equation following the original idea due to Hinze. As an advan- 
tage, we can immediately isolate the contributions due to the mean shear to identify the 
range of scales where they become the leading terms. Next, by generalizing the concept 
of structure functions, we are able to study systematically the scaling properties of shear 
turbulence and its universality with respect to homogeneous and isotropic conditions. 

The Kolmogorov equation for homogeneous isotropic turbulence follows from the Karman- 
Howarth equation 

< 5V? >= --er + 6iy-^ < SV? > (1.1) 

" 5 dr " 

where SV\\ = [u{x-\-f) — u{x)] -f/r are the longitudinal velocity increments. In the inertial 
range, the Karman-Howarth equation reduces to the well known four-fifths law relating 
the third order structure function to the mean rate of energy dissipation e and to the 
separation r. In absence of intermittency effects the four-fifths law suggests that, beyond 
the separation r, any inertial range scalings should involve only the mean rate of energy 
dissipation e, or, more properly, the average energy flux through the inertial range. In 
other words, scalings in terms of separation for structure functions of order p should 
necessarily be given by 

< dV( >oc eP/\P/^ . (1.2) 

It follows that the dimensionless ratios of structure functions must go to a constant as 
the separation is reduced within the inertial range 

lim kSV,? >/< SV? >P/^ ^ const , (1.3) 
i.e. no intermittency is present. This is essentially the K41 theory of homogeneous 



isotropic turbulence (see Frisch). Preliminarily we may observe that eq. (1.3) implies 
< SVr >oc< SV? >p/i. 



When intermittency is accounted for, eq. (1.3) is not valid. As a matter of fact scaling 



laws are nonetheless observed in the inertial range, with exponents substantially different 



from the dimensional prediction for large p (Benzi et al. 199£). The observed behavior 
can be recast in the form 



< 5V[ >oc j^-^ j , (1.4) 

where Lq is the integral scale of turbulence and t(p/3) are the intermittency corrections 
to the dimensional scaling exponents. 



Jhi the context of the Kolmogorov Refined Similarity Hypothesis (RKSH) ( Kolmogorov 



1962| ) the intermittency corrections are associated to the statistical properties of the 
coarse grained energy dissipation, 

IB. 



Br 



Qoc dVr, (1.5) 
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with eioc the local rate of energy dissipation, by expressing the structure functions as 

< >cx< eP/^ > rP/\ (1.6) 

In words, the intermittency of tlie velocity field is originated by an intermittent energy 
cascade which ends up in an intermittent dissipation field. The link between the two 
fields is provided by eq. ( |1.6D where 



(1.7) 



Intermittency effects on vel ocity structure fu nctions in the inertial range are documented 
beyond any doubt (see e.g. Zhou fc Antonia ). More controversial is the assessment of eq. 
(1.6). However the evidence is all in favor of tliis equation, which, though not derived 
from first principles, seems to describe correctly the physics of small scale intermittency 
( Wang et ai ). In this context, the four- fifths law is an exact constraint to be fulfilled 
by any acceptable theory of turbulence. Besides this, the Kolmogorov and the Karman- 
Howarth equations may also suggest the basic ingredients of a reasonable turbulence 
theory. 

We can now re-express more technically the aim of the present paper. Namely, starting 
from the extension of the Karman-Howarth equation to the homogeneous shear flow (see 
also pbcrlack for a related derivation), we are interested in generalizing the RKSH to 
anisotropic turbulence. As we shall see, the proposed similarity law includes the classical 
Kolmogorov- Obhukhov theory of intermittency as a special case. Further, we intend to 
check the proposed form of similarity by using DNS results of a homogeneous shear flow. 
This may make things appear a little less straightforward, since, given the limitations 
of DNS, scaling laws cannot be detected directly in terms of separation. To overcome 
the difficulty we analyze our data via Extended Self Similarity ( Benzi ct al. 199^ ), to 
show that, indeed, the generalized RKSH holds for the homogeneous shear flow. Actually 
ESS scalings, < dVf >oc< SV.-' >^'-p\'^^ with (:{p\q) / p/q, are observed also at moderate 



Reynolds number both for isotropic and for shear turbulence (see Antonia et al, Ruiz 
Chavarria et al. and Toschi ct al). In particular, for homogeneous isotropic turbulence 



and within the ESS rationale one can rewrite eq. (1.6) as 



< 6VP >cx < >/eP/^ < SVf >p/^ , 



(1.8) 



to check its validity against moderate Reynolds number data (Benzi et al. 1996). 

The generalization of the Kolmogorov-Obhukhov theory proposed here is well sup- 
ported by the data and rationalizes a number of recent results concerning intermittency 
in shear dominated turbulence. Actually a number of papers have recently appeared on 
this subject. Strictly related to the present results is the paper by Benzi et al. (1999) 
where a new form of similarity law has been proposed to explain the increased inter- 
mittency of the near wall region of wall bounded flows noticed by a number of authors 
(see Ruiz Chavarria et al., Toschi et al.). That form of similarity law has been checked 



by using experimental (Jacob et al.) and DNS data (Gualtieri et al.), to finally address 
the possible cohexistence of two different intermittent regimes (Casciola et al.). In the 
present paper we discuss a unifying formulation showing that the classical RKSH and 



the form proposed by Benzi ct al. (1999) are obtained as limiting cases of a more general 
similarity law. 

Despite the fact that the Karman-Howarth equation for homogeneous isotropic turbu- 
lence is conveniently expressed in terms of third and second order longitudinal structure 



4 



C. M. Casciola, P. Gualtieri, R. Benzi and R. Piva 



functions, more generally it involves integrals on the sphere of radius r such as 

< SV^SV\\ > dSr (1.9) 
etc. The primary form of Karman-Howarth equation is recast in the more familiar form 



(1.1) by using isotropy. Clearly, the last step can not be applied to anisotropic turbulence. 
Consistently the generalized Karman-Howarth equation gives a relation between integrals 
on the sphere, implying that the generalized RKSH will also maintain the same flavor. 



2. Karman-Howarth equation 

In terms of fluctuating components, u, a homogeneous field v may be written as 

dUi 

Vi=Ui + —-^Xk. (2.1) 

OXk 

Here the velocity field obeys the incompressible Navier-Stokes equations and the mean 
component is solenoidal, dUi/dxi = 0. Standard notations are used for quantities evalu- 
ated at X and a prime denotes those a,t y = x + r. After multiplying the equation for Ui 
by u'j and that for u'^ by Ui, adding the results and averaging we obtain the equation for 
the two points correlation tensor Rij{f,t) =< Uiu'^ >. The equation for the trace 
reads 

-^{R.k+Rk.) + 2i.p^, (2.2) 
OXk orkOTk 

where, by homogeneity, d/dxk < • >— —d/drk < • > and d/dyu < • >= d/drk < • >, 
and the pressure term drops due to incompressibility. Equation ( p.2[ ) can be re-expressed 
in terms of structure functions, by considering the velocity increments 

SVt{x, r) = Ui{x + f) - Ui{x) 

(2.3) 

5V\\{x,r)^5V,-, 
r 

SV^ = SVtSVi, and 5Uk = 11'^. -Uk- Actually, since R.^i =< UiUi > -5 < SV'^ > and 
Iii,k + Rk,i UiUi > — < 6ViSVk >, and considering the standard manipulations 

< SV'^SVk >= — 2[< Uiu'^u'k > — < Uiu[uk >]+ < UiUiu'f. > — < u^u'^Uk > , (2.4) 

the Karman-Howarth equation for the velocity increments follows as 

d dU- dU d"^ 

— < 6V^SVk+SV^6Uk > +2—^ < SV^SVk 4—^ < u,Uk > +2i^— — < SV^ >, 
drfc dxk OXk orkOTk 

(2-5) 

where use has been made of incompressibility. Here the first term in the right hand side 
is identified with the average dissipation rate, given the balance between production of 
turbulent kinetic energy and viscous dissipation, 

dUi dui dui 

- — < UiUk >= -V < - — - — >= -e. (2.6) 

OXk OXk OXk 



All averaged terms in eq. (2.5) are only function of separation and integration over the 
sphere of radius r yields 



< SV^SVw > + < 5V'^5U\\ > dSr 



dB,. 



8U- 

2S^ < SV^SVk > dVr 



= -4er - 



dr 



dBr 



2v < SV^ > dSr . (2.7) 



When U = {Sy, 0, 0), by denoting by ni and n2 the components of the normal in the 
directions of mean flow and mean gradient, respectively, we finally arrive at 



1 



47rr^ 



< SV'^SVn > +nin2Sr < 6V^ > dS,. 



dBr 



2S < 6u6v > dVr 



4_ 1 d 

€T -\- 

3 47rr2 dr Jgg^ 



2v < SV^ > dSr . (2.. 



3. Scale by scale budget 

In homogeneous isotropic turbulence the Karman-Howarth equation describes the pro- 
cess of energy cascade towards small scales and establishes the balance between energy 



flux and dissipation according to the four-fifth law (l.I). The energy, injected at the in- 



tegral scale Lo, is simply transfered across the inertia! range to be dissipated at scales of 



the order of the Kolmogorov length r]. Equation (2^) extends this result to anisotropic 
flows where the cascade process is strongly modified by the continuous energy injection 
associated to the mean velocity gradient. 

Following the derivation of the previous section, the non linear terms generate two 
distinct contributions. In particular the correlation 

< SVtl >= — ^ / < SV'^SVii > +nin2Sr < SV^ > dSr (3.1) 

where the subscript tr stays for transfer, is the proper generalization of the third order 
structure function < SV^^ > to non isotropic flows. Its meaning is rather clear: the flux 
across the surface dBr arise due to transport of turbulent kinetic energy operated by 
both the mean flow Sr and by the turbulent fluctuations SV^^ . It arises from the terms in 
divergence form in eq. ( |2.5| ), as appropriate for a flux. The second term in the left hand 
side of eq. ( p.Sf ), namely 

< ^^pr >= TZ3 / 25* < 6uSv > dVr, (3.2) 



Br 



where the subscript pr stays for production, describes the rate of turbulent kinetic energy 
production up to the scale r and keeps its form as a volume integral. 

All the terms appearing in eq. (2.8) are plotted in figure |l| and clearly sum up to zero 



for all separations r. In the right part of the figure, transfer < SV^^. > and production 
< SVpr > are compared to — 4/3er deprived of the finite Reynolds number contribu- 
tions. At the large scales e is balanced by production. At small scales, instead, e is 
balanced by energy transfer. Actually, moving towards the scales of the classical inertial 
sub-range from above, the production having place in the larger scales piles up and orig- 
inates the flux of energy which, transfered down the cascade, is eventually dissipated. In 
fact the production and the transfer terms achieve the same order of magnitude at the 
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Figure 1. Left: The different terms of eq. (2^) are plotted against r/rj. < SVtr >, triangles. 



< SVpr >, circles, the viscous correction, squares, and — 4/3er diamonds. The four terms sum 
to zero, solid line. Right: The squares denote — 4/3er + d/dr Iv < &V'^ > dSr, the other 
symbols are defined as for the left part of the figure. 



characteristic length scale where < (JV^"^ >~< W^^ >. The cross over scale may be esti- 



er and < SV^^ >c 



SrSuSv with 



mated on dimensional grounds by assuming < SVf^ >: 
Su 0:^ Sv = e^/^r^/^^ to yield the classical expression of the shear scale Ls — y^e/S^. This 
is actually well supported by the data shown in figure |l| where the balance is achieved at 
rb/rj ~ 25, roughly corresponding to the dimensional prediction of Lg/ri c:± 20. 



4. Similarity laws 

The generalized Karman-Howarth equation (2.8) described in the previous sections, in 
the limit of large Reynolds numbers, reduces to 



< SV^ >oc er. 



(4.1) 



where 

6V^^6V,l + 6V;^,. (4.2) 

This equation is the analogous for the shear flow of the four-fifths law of homogeneous 
isotropic turbulence, to which it nicely reduces for zero shear. On purely dimensional 
grounds, it should generalize to higher order moments as < SVp >(x eP/^r-f/^. To incor- 
pora te in termittency corrections we follow the line of reasoning th at le d to the similarity 
law (L6) by considering the coarse-grained dissipation field, eq. (1^), as the carrier of 
the intermittency corrections to the dimensional scalings of structure functions. In our 
case this procedure leads to the prediction 

< SVP >oc< > (4 3) 

Unfortunately we cannot validate equation ( |4.3| ) directly, since, given the limitations 
in the Reynolds number of our DNS, no power-law behavior can be observed when using 
the separation r as scaling variable. However, the concept of Extended Self-Similarity 
can help overcoming the difficulty by recasting our ansatz in the form 



(4.4) 
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Figure 2. Left: < SV^ > vs. < X 6V^ >^ (open circles). Data are fitted in the whole range 
of scales by a power law with slope s = 1 (solid line). The shear scale approximately corresponds 

to the sixth point from the left. Right: The ratio < SV^ > /(< >< SV'^ >^)" vs. separation 
r/rj for a — 1.05 (triangles), a = 1.00 (circles), a = 0.95 (squares). 



Whenever equation (4.3) holds, equation (4.4) follows by using (4.1) to eliminate r in 
favor of < SV^ >, though ESS should work under more general conditions than ( [4.3[). 

As shown in figure ^ for p = 6, the data we have available well support equation (4.4) 
which apparently holds uniformly through the scales. This is a peculiar feature of the 
present description, where all aspects of shear turbulence are retained, both the produc- 
tion term associated to SVp^. and the transfer term associated to SV^j. . The two contribu- 
tions are evaluated exactly avoiding any approximation. 

As a further check, we have also reported in figure ^ the ratio < SV^ > /(< >< 
SV^ >2)" for different a's. The value a = 1, giving eq. (4.4), achieves the optimal 
compensation around a constant for the whole range of scales. Quantitatively, the com- 
pensated data for eq. (4.4) fluctuate about 4.5% giving clear evidence in favor of the 
proposed similarity law above and below the shear scale Ls- 

According to eq. (4^), the random function 6V^ is statistically equivalent to the coarse 
grained energy dissipation e^. It is a fact that, for the present data where the shear affects 
the turbulent fluctuations enough to break the classical refined similarity, the statistics 
of the dissipation fleld shows no appreciable difference with homogeneous isotropic tur- 
bulence. In figure ^ we compare the second order moment of for the homogeneous 
shear flow (solid line) and for homogeneous isotropic turbulence (triangles). The corre- 
spondence between the two statistics is confirmed by analyzing higher order moments, 
such as the third order one plotte d o n the right of the same figure. 

The above observation and eq. ( 4^ ) allow for the conclusion that the generalized struc- 
ture functions < 5V > are unaffected by the shear, see the left part of figure^. Actually, 
the relative scaling exponents of < SVp > vs. < SV^ > are, within 1% accuracy, the same 
found in isotropic conditions. Under this respect, they extend to shear dominated fiows 
the universal behavior well documented in homogeneous isotropic turbulence. 

As an additional comment, we like to compare the present findings with related work 
concerning the similarity law proposed by Benzi et al. (1999)] for shear dominated flows 
and the double intermittent behavior in shear turbulence discussed by Casciola et al . For 
r < Lg, the budget provided by the Karman-Howarth equation suggests that SV^"^ :$> 5V^^. 
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Figure 3. Left: < > vs. r/77 in homogeneous isotropic turbulence (triangles) and in the ho- 
mogeneous shear flow (solid line). Right: < > vs. < > in homogeneous isotropic turbulence 
(dashed line, slope s = 2.92) and in the homogeneous shear flow (open circles). 



In this limit eq. (4.4) reduces to the classical Kolmogorov-Obhukhov similarity, 

< SVl >oc< >< SVtl >P'^ . (4.5) 
This equation is plotted for p = 6 on the right of figure 0, where it is shown to hold 



below the shear scale as indicated by the solid line. Equation (4^) still involves integrals 
over the sphere. It may be recast into the familiar form involving only the longitudinal 
increments, evaluated in the streamwise direction 5u. This is done by assuming that, as in 
the homogeneous isotropic case, 5Vtr — A{r)5u where A{r) is a suitable non fluctuating 
function. Hence eq. (4.5) is re-written as 

< 5uP >oc< >< 5u^ >P/^ . (4.6) 

On the contrary in the range of scales r > Lg the scale-by-scale budget shows that the 



only productions terms in cq. (4.4) are relevant. In this range we obtain asymptotically 

>P/3 (4J) 



>oc< eP/^ >< SVj^^ 



Equation (4.7) is also plotted for p = 6 on the right of figure ^ and manifestly holds at the 
largest scales. Equation ( 4.7) exactly corresponds to the similarity law originally proposed 
by Benzi ct al. (1999). Actually, on dimensional grounds, 5Vp^ may be approximated 
in terms of longitudinal and transverse velocity increments in the streamwise direction 
SVpj, ~ B[r)5u5v where B{r) is, again, a non fluctuating function. Hence, by assuming 
that 5u5v ~ 5v? , the similarity law 



< 5vF >(x< >< 5u^ >f/2 



(4.8) 



follows. As discussed in the introduction, equation (4.S) has been repeatedly verified, both 
numerically and experimentally, in a number of different shear flows and has been shown 
able to correctly capture the change in the nature of the intermittency of longitudinal 
increments across the shear scale. 
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Figure 4. Left: < SV^ > vs. < SV^ > in homogeneous isotropic turbulence (dashed 
hne with slope s — 1.78) and in the homogeneous shear flow (open circles). Right: 
< 5V^ > vs. < >< 5V^ Production (a = pr) filled circles, transfer {a — tr) open 
circles. The slope of the two solid lines is s = 1. 



5. Final remarks 

The exact Karman-Howarth equation for the homogeneous shear flow has been ex- 
ploited to provide a rational framework for the analysis of intermittency in shear tur- 
bulence, i.e. in flows characterized by large scale anisotropics. Besides conflrming that 
production of turbulent kinetic energy prevails, as expected, at large scales while transfer 
is crucial at the smallest ones and that the two ranges are separated by the shear scale 
Ls, the approach has proven fruitful in allowing the generalization of the Kolmogorov- 
Obhukhov theory of intermittency to shear dominated turbulence. In particular the two 
main achievements of the paper are summarized as follows. A generalization of the struc- 
ture functions can be introduced in shear turbulence, to extend the classical description 
in terms of velocity increments used in the standard theory for isotropic flows. The gen- 
eralized structure functions obey a similarity law which, formally, keeps the structure of 
the Kolmogorov-Obhukhov similarity, in such a way that in the limit of zero shear the 
classical results of isotropic turbulence are recovered. This generalized similarity law has 
been favorably checked against DNS data of a homogeneous shear flow and it has been 
shown to be consistent with both the classical reflned similarity at small scale a nd the 
new form of similarity law recently introduced for shear dominated turbulence (Benzi 
et al. (1999) ). While scaling laws in terms of traditional longitudina l structure func tions 



display two different forms of intermittency across the shear scale (Casciola et ai), the 



generalized similarity described in this paper holds uniformly across the entire range of 
available scales. 

As a further contribution, the properties of the dissipation field related to intermit- 
tency have been shown to be unaltered by the shear, i.e. they are identical to those of 
isotropic turbulence. This implies that the intermittency of the dissipation field displays 
universality, which entails universality of the generalized velocity increments we have de- 
fined. By exploiting the relationship between the generalized and the standard velocity 
increments we find that the universal intermittency corrections for the generalized incre- 
ments give reason for the increased intermittency of the standard velocity increments, 
in agreement with previous numerical and experimental results in shear turbulence. In 
particular, since the longitudinal structure functions of order p correspond to generalized 
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moments of the order p and 3/2 p below and above the shear scale, respectively, the flat- 
ness factor for the standard increments increases across the shear scale as a consequence 
of the universality manifested by the generalized structure functions. 

In conclusion we like to point the attention of the reader on the potential interest of 
the present work, providing a uniform description for both shear dominated and isotropic 
turbulence, in the context of LES of wall bounded flows. Actually the present descrip- 
tion provides an appropriate estimate for the energy flux in non isotropic conditions, a 
crucial point in developing working closure in LES ever since the now classical model of 
Smagorinski. It could be also worth emphasizing the role of the homogeneous shear flow 
as the proper setting to check possible alternative formulations of LES closures. 
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